Quasi-one-dimensional anisotropic Heisenberg model in a transverse magnetic field 
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The phase diagram of weakly coupled XX Z chains in a transverse magnetic field is studied using 
the mean-field approximation for the interchain coupling and known exact results for an effective 
one-dimensional model. Results are applied to the quasi-one-dimensional antiferromagnet CS2C0CI4, 
and the value of interchain interaction in this compound is estimated. 



The effects induced by magnetic fields in low- 
dimensional magnets are subjects of intensive theoretical 
and experimental research 1]. One of the striking effects 
is the dependence of magnetic properties of quasi-one- 
dimensional (QID) antiferromagnets with anisotropic in- 
teractions on the direction of the applied magnetic field. 
For example, the behavior of these systems in a trans- 
verse magnetic field is drastically different in comparison 
with the case of the longitudinal field applied along the 
anisotropy axis. In particular, the transverse field in- 
duces a gap in the spectrum and the antiferromagnetic 
long range order (AF LRO) in the perpendicular direc- 
tion. A quantum phase transition takes place at some 
critical field, where the LRO and the gap vanish. The 
phase transition of this type has been observed in the 
QID antiferromagnet CS2C0CI4 The simplest model 
of the one-dimensional anisotropic antiferromagnet in the 
transverse field is the spin-i XX Z chain described by the 
Hamiltonian 

Hid = J ^{SnSn+i + S^S^^^ + AS^S^^^) ~ H^S^^ 

(1) 

where A is an anisotropy parameter which assumed to 
be < A < 1. 

It was proposed that low-energy properties of 
CS2C0CI4 is described by with J ~ Q.23mev and 
A = 0.25. In contrast to the case of the longitudinal field 
the symmetry-breaking transverse field does not com- 
mute with the XX Z Hamiltonian and the exact integra- 
bility of 1^ is destroyed. The model has been investi- 
gated using diS^erent approximate approaches ^ H, IS - 
The scaling estimates at small field show that the 
transverse field generates the staggered magnetization 
M^t = {{-!)" Sy) (AF LRO in the Y direction) and the 
gap in the spectrum m (at H = the spectrum is gapless) 
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To study the model (Q, when the field H is not small, 
the mean-field approximation (MFA) has been proposed 
in and elaborated in 0. The MFA is based on the 
Jordan- Wigner transformation of spin-i operators to the 
Fermi ones with the subsequent mean-field treatment of 
the four-fermion interaction term. As a result the aris- 
ing Hamiltonian is quadratic in Fermi-operators and it 
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FIG. 1: The dependence of LRO parameter on magnetic field 
for ID chain (dashed line) and quasi-lD system (solid line) 
for A = 0.25. 



is solved exactly. Transforming this MFA Hamiltonian 
back to spin variables we obtain a spin-i XY model in 
the longitudinal field 



(3) 

where parameters J', 7 and h are determined by the MFA 
self-consistent conditions 

The model © is exactly solvable and its properties 
are well studied 9J. This model undergoes the T = 
phase transition of the 2D Ising universality class at 
h = \ corresponding to the MFA value of the critical field 
i/j!^^(A). In particular, in the vicinity of the critical field 

Mst ^ \Ul^{,A) - A comparison of the MFA re- 

sults with those obtained in precise numerical DMRG cal- 
culations shows high accuracy of the MFA at ^ J |3| ■ 
The dependence Mst{H) for A = 0.25 obtained with use 
of the MFA and scaling estimate (0) is shown on FigQJby 
dashed line. This magnetization curve is qualitative sim- 
ilar to that observed in neutron-scattering experiments 
on CS2C0CI4. At the same time, there is essential differ- 
ence in the low-field behavior of Mgt- The experimental 
AF ordered moment is finite at H = while Mst ac- 



2 



r+S 



ay qV 
^n,r^7i+l,r 

ay qV 



(4) 



cording to This difference is due to weak interchain 
couphngs in real systems and these couphngs form a 3D 
magnetically ordered moment below a Neel temperature 
T/v- Besides, interchain couplings extend the ID ordered 
phase with Mgt ^ to finite temperatures. Therefore, 
to describe low temperature properties of real QID com- 
pounds it is necessary to take into account interchain 
interactions. 

In this latter we will consider the system of coupled 
parallel XX Z chains in the transverse field described by 
the Hamiltonian 

-J.I 

where n and r label lattice sites along the chain and in 
perpendicular directions, S is summed over two nearest 
neighbor vectors in the transverse directions, Jx is a weak 
coupling between neighboring chains. 

A standard method for treating the model l@J is to use 
the mean-field approximation for interchain coupling and 
to treat the resulting effective ID problem as exactly as 
possible [13, EH (we call this approach as chain mean- 
field theory (CMFT) to distinguish it from the MFA for 
the ID model (^)). We assume that AF order in each 
chain to be oriented along the Y direction and the uni- 
form magnetization along the X axis as it occurs in the 
pure ID model The quasi-lD model contains another 
mechanism to generate the LRO. If one of the chains is 
AF ordered, the interchain couplings induce an effective 
staggered field on the nearest chains. In the CMFT in- 
terchain coupling is replaced by effective fields and the 
Hamiltonian Q reduces to an effective ID Hamiltonian 
having the form 

-{H-H^)Y,Sl^HyY,{-irSl (5) 

where fields and Hy are determined by self- 

consistency relations 



(6) 



z is the transverse coordination number. 

At first, we consider the model Q at _ff = and T — Q. 
It can be easily shown that the self-consistency relation 
gives (S'^) = and the model © reduces to the XXZ 
chain in the staggered field. The low-energy properties 
of this model are described by a quantum sine-Gordon 
model [13 

n = Ho + V 

-Ho = dx{{d,<i>f + (9,9)2} 

V = -Hy^y2A{Tj) j dxsi^{y/2^Q) (7) 



where $(a;) and Q{x) are boson and dual fields respec- 

Jsin{7rr;) 

efficient ^(77) was found in llj. 



tively, v{vi) = 2-2l^' the sound velocity and the co- 
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The spectrum of Ti.Q is gapless. The perturbation V 
has scaling dimension 77/2 and generates the mass gap 
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where constant C is \IM 
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The sta gge red magnetization Mst is related to a mass 
gap m as [13 
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From the equations and we get 



M,, = D[CD 



T7/2 



m = v\CD 



zJ± 



(13) 



(14) 



The AF LRO M^t survives at T < Tn- The Neel 
temperature T^r can be found using the random phase 
approximation (RPA) . The RPA dynamical susceptibility 
of coupled chains in disordered phase (T > T^) is 



X^nc.,fc,fci) = 



The condition determined Tjv is 



(15) 



zJxX3'^(0,^) = l (16) 

The dynamical susceptibility of the ID XXZ model 
at T < J is known 113 



B / V \2-r; 



xS(0,») = -(5^) 



(17) 
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where 

Using the condition H16() we extract the Neel tempera- 
ture at iJ = 



2n \ V 



(19) 



0.0147 



We note that the ratio Tn/ui does not depend on J± 
and is determined by ID parameter rj only. 

An analysis of experimental data carried out in Rcf. 2] 
has shown that the quasi-lD antiferromagnet CS2C0CI4 
consists of two interpenetrating sublattices with iden- 
tical intrasublattice interactions. These sublattices are 
non-interacting on the CMFT level. Each sublattice 
has tetragonal symmetry and described by the model 
Q with z = 4. However, no direct experimental data 
on the value of the interchain interaction Jx is avail- 
able. The Neel temperature in CS2C0CI4 at iJ = is 
Tn = 0.0813 J = 0.217i<: 0. Using these data we can es- 
timate unknown value of J± in CS2C0CI4. Substituting 
A = 0.25 {A = 0.1405) in (O we find 

J± 
J 

This value is really small, so that our assumption about 
quasi-lD behavior of the system is justified. 

Further, using the found value of J± we can find the 
staggered magnetization Mst at T = 0. According to 
Eq.lO Mst = 0.348. The experimental value of the AF 
ordered moment at T < Ta? is Mst ~ 0.342 2]. Such a 
perfect coincidence confirms our estimate (|20|l . Besides, 
the found value of J± gives us also the gap (tn|) m — 
0.78K. It is remarkable that even so small interchain 
coupling as in Ea. l|^ causes so large value of LRO. 

At = and T = the AF LRO is generated by 
the interchain couplings. At _ff > the 'one-dimensional' 
mechanism is switched. The crude estimation of the value 
if*, at which this mechanism becomes predominant, can 
be obtained by a comparison of (0) with H14|l 



(20) 



H* 



JizJ^/J) — 



(21) 



At H > H* in the Hamiltonian © the mean field Hx 
can be neglected in comparison with H and at T = the 
main effect of Hy consists in a small shift of the critical 
field H^^ (see below). 

At H = and Hy = the spectrum of the 

model © is gapless. The perturbation HyJ2{—^)^S^ 
has scaling dimension 1/8 and generates the mass gap 
TO ~ {Hy/Jf/^^ and AF LRO M^t ~ {Hy/J)^/^^ in the 
model (dJ. The self-consistency relations ^ therefore 
give 



1/14 



Mst{H^n ~ izJ±/J) 



(22) 



To estimate the Neel temperature Tm{H) in the RPA 
it is necessary to know the finite temperature staggered 
susceptibility Xioi^' ^'^^ model ^ at H > 0. Un- 
fortunately, it is unknown. Instead, we consider the MFA 
model ©, for which the susceptibility can be found. As it 
was noted above the MFA describes correctly the ground 
state properties of the model ^ at H > J. We expect 
that the MFA gives a satisfactory description of at 
low temperature (T <C J) as well. The problem of find- 
ing Tm{H) can be solved in the same manner as it was 
done by Carr and Tsvelik in |l6| for quasi -ID quantum 
Ising model, which on the CMFT level reduces to (O 
with 7 = 1. 

We are mainly interested in the region of the fields 
near the critical field H^^{A) or at ft- ^ 1 in terms of 
the MFA Hamiltonian Exactly at h = 1, where the 
model © is critical, the staggered susceptibility at T ^ 
J according to 15] is 



^^^KS^^t) F (1/8) r2 (15A6) ^^^^ 



where sound velocity at the critical field — 7J' is de- 
termined from the MFA self-consistent equations and 
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with Glaisher constant A ~ 1.282. 
The Neel temperature in the RPA is 
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For A = 0.25 the critical field in the MFA is Hl'^ w 
1.6 J and Vc ~ 0.185 J. Therefore, the estimated Neel 
temperature for Cs2CoCh is Tn{HI°) = 0.145i4:. 

Near the critical field at > H],^ (disorder region in 
ID model) the low-temperature staggered susceptibility 
is well approximated by the formula 1] 
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where the gap ra — \H ^ 
time r, 

phase transition (|16|l reads 



1 + 7 v'^k^ + m'^ ~ {uj + i/tc) 

' and the phase relaxation 



^^m/T ^ In this case the RPA condition of the 
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At first we estimate the shift of the critical field SH^ — 
H^^ — caused by interchain couplings. This shift is 
determined by the condition T ^ in Eq. 127|l 
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FIG. 2: Expected phase diagram of the model 

For A = 0.25 and found value of Jj. ((201) the shift of the 
critical field is about 3%. We note that in the vicinity of 
the critical point H^'^ the low-energy properties of QID 
model Q belong to the universality class of the (3+1)- 
dimensional classical Ising model. 

Eq.|j23) gives also the behavior of Neel temperature 
near the 3D critical point H^^ — iJ <C 5Hc 

T^{H) « 25H, In-i (^ Joljj ) (29) 

At intermediate fields H* < H < H}^ (ID ordered 
region) the low temperature T <^ m staggered suscepti- 
bility has an exponential form Q 

„3/4 

xZiO, ^) - {mlvfl^T, ~ ^1/4^3/2 ^""^"^ (30) 

with correlation length = "x/i^e™^^ 'I'l- Thus, for 
zJ^ <C m{m/vY/'^ the RPA criteria yields 

Tm{H) ^ 2m{H) In-i (31) 

Combining the found expressions for Neel temperature 
in different regions (^HJ , , 123) , IpH^ we arrive at the 
phase diagram schematically shown on FigI21 Since the 
gap m{H) in the AF ordered region has a maximum at 
some intermediate value of field .7^^ 8], then according 
to Ea. l|31() the function Tm{H) also has a maximum as 
shown on Fig|H This fact was experimentally observed 
in Cs2CoCIa ^\. 

In summary, the main effect of interchain coupling on 
the system of coupled XX Z chains consists in extending 



of the AF ordered phase to finite temperatures. Based 
on known exact results for the effective QID model we 
estimated the value of interchain coupling in CS2C0CI4, 
which turns out to be very small. We found that even 
weak coupling between XX Z chains can cause a remark- 
ably large AF LRO below the Neel temperature. 

In the presence of transverse magnetic field the 3D ef- 
fects become less pronounced. They lead to a small shift 
of the critical field and to a small correction to AF LRO. 
The found small interchain coupling in CS2C0CI4 im- 
plies that the properties of this compound in transverse 
magnetic field are similar to that of ID quantum Ising 
model in a magnetic field. The spectrum of this model 
at exactly solvable critical point consists of eight particles 
which may be observed as a few coherent peaks in the dy- 
namical magnetic susceptibility of Cs2CoCli in the AF 
region slightly below the observed critical field [l^ . 

This work was supported under RFBR Grant No 03- 
03-32141. 



* Electronic address: |dmitriev@deom.chp h.ras.ru| 
[1] S.Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, UK, 1999). 
[2] M.Kenzelmann, R.Coldea, D.A.Tennant, D.Visser, 
M.Hofmann, P.Smeibidl and Z.Tylczynski, Phys. Rev. B 
65, 144432 (2002). 
[3] J.Kurmann, H.Tomas and G.MuUer, Physica A112, 235 

(1982) ; G.MuUer and R.E.Shrock, Phys. Rev. B32, 5845 
(1985). 

[4] S.Mori, J.-J.Kim and I.Harada, J. Phys. Soc. Jpn 64, 
3409 (1995). 

[5] Y.Hieida, K.Okunishi and Y.Akutsu, Phys. Rev. B64, 

224422 (2001). 
[6] A.Langari, Phys. Rev. B69, 100402 (2004). 
[7] D.V.Dmitriev, V.Ya.Krivnov, A.A.Ovchinnikov and 

A.Langari, JETP 95, 538 (2002); D.V.Dmitriev, 

V.Ya.Krivnov and A.A.Ovchinnikov, Phys. Rev. B65, 

172409 (2002). 

[8] J.-S.Caux, F.H.L.Essler and U.Low, Phys. Rev. B68, 
134431 (2003). 

[9] E.Barouch and B.M.McCoy, Phys. Rev. A3, 786 (1971). 
[10] D.J.Scalapino, Y.Imry and P.Pincus, Phys. Rev. Bll, 
2042 (1975). 

[11] H.J.Schulz, Phys. Rev. Lett. 77, 2790 (1996). 
[12] I.Affleck and M.Oshikawa, Phys. Rev. B60, 1038 (1999). 
[13] S.Lukyanov and A.Zamolodchikov, Nucl. Phys. B493, 
571 (1997). 

[14] Al.B.Zamolodchikov, Int. J. Mod. Phys AlO, 1125 (1995). 
[15] H.J.Schulz and C.Bourbonnais, Phys. Rev. B27, 5856 

(1983) . 

[16] S.T.Carr and A.M.Tsvelik, Phys. Rev. Lett. 90, 177206 
(2003). 



